#include<bits/stdc++.h>
using namespace std;
namespace IO{
    const int BUF=1<<21;
    char buff[BUF],*p1=buff,*p2=buff;
    #define gc() (p1==p2&&(p2=((p1=buff)+fread(buff,1,BUF,stdin)),p1==p2)?EOF:*p1++)
    template<typename T>inline void read(T &x){
        char ch=gc();x=0;int f=1;
        while(!isdigit(ch)){if(ch=='-')f=-1;ch=gc();}
        while(isdigit(ch)) x=(x<<1)+(x<<3)+(ch^48),ch=gc();
        x*=f;
    }
}
using namespace IO;
const int N=3005;
int n,m,l[N][N],r[N][N],u[N][N],d[N][N],tag[N],cnt1[N][N],cnt2[N][N];
long long ans;
char s[N][N];
void solve(int lt,int rt,int L,int R){
	if(lt>=rt||L>=R) return;
	if(rt-lt>R-L){
		int mid=(lt+rt)/2;
		for(int i=L;i<=R;++i){
			for(int j=max(u[mid][i],lt);j<=mid;++j) ++tag[min(r[j][i],R)];
			int res=0;
			for(int j=R;j>=i;--j){
				res+=tag[j],tag[j]=0;
				if(u[mid][i]>=u[mid][j]) cnt1[i][j]=res;
			}
		}
		for(int i=R;i>=L;--i){
			for(int j=max(u[mid][i],lt);j<=mid;++j) ++tag[max(l[j][i],L)];
			int res=0;
			for(int j=L;j<=i;++j){
				res+=tag[j],tag[j]=0;
				if(u[mid][i]>u[mid][j]) cnt1[j][i]=res;
			}
		}
		for(int i=L;i<=R;++i){
			for(int j=min(d[mid+1][i],rt);j>mid;--j) ++tag[min(r[j][i],R)];
			int res=0;
			for(int j=R;j>=i;--j){
				res+=tag[j],tag[j]=0;
				if(d[mid+1][i]<=d[mid+1][j]) cnt2[i][j]=res;
			}
		}
		for(int i=R;i>=L;--i){
			for(int j=min(d[mid+1][i],rt);j>mid;--j) ++tag[max(l[j][i],L)];
			int res=0;
			for(int j=L;j<=i;++j){
				res+=tag[j],tag[j]=0;
				if(d[mid+1][i]<d[mid+1][j]) cnt2[j][i]=res;
			}
		}
		for(int i=L;i<=R;++i)
		for(int j=i+1;j<=R;++j){
			if(s[mid][i]==s[mid+1][i]) ans+=1ll*cnt1[i][j]*cnt2[i][j];
			cnt1[i][j]=cnt2[i][j]=0;
		}
		solve(lt,mid,L,R),solve(mid+1,rt,L,R);
	}else{
		int mid=(L+R)/2;
		for(int i=lt;i<=rt;++i){
			for(int j=max(l[i][mid],L);j<=mid;++j){
				++tag[min(d[i][j],rt)];
			}
			int res=0;
			for(int j=rt;j>=i;--j){
				res+=tag[j],tag[j]=0;
				if(l[i][mid]>=l[j][mid]) cnt1[i][j]=res;
			}
		}
		for(int i=rt;i>=lt;--i){
			for(int j=max(l[i][mid],L);j<=mid;++j) ++tag[max(u[i][j],lt)];
			int res=0;
			for(int j=lt;j<=i;++j){
				res+=tag[j],tag[j]=0;
				if(l[i][mid]>l[j][mid]) cnt1[j][i]=res;
			}
		}
		for(int i=lt;i<=rt;++i){
			for(int j=min(r[i][mid+1],R);j>mid;--j) ++tag[min(d[i][j],rt)];
			int res=0;
			for(int j=rt;j>=i;--j){
				res+=tag[j],tag[j]=0;
				if(r[i][mid+1]<=r[j][mid+1]) cnt2[i][j]=res;
			}
		}
		for(int i=rt;i>=lt;--i){
			for(int j=min(r[i][mid+1],R);j>mid;--j) ++tag[max(u[i][j],lt)];
			int res=0;
			for(int j=lt;j<=i;++j){
				res+=tag[j],tag[j]=0;
				if(r[i][mid+1]<r[j][mid+1]) cnt2[j][i]=res;
			}
		}
		for(int i=lt;i<=rt;++i)
		for(int j=i+1;j<=rt;++j){
			if(s[i][mid]==s[i][mid+1]) ans+=1ll*cnt1[i][j]*cnt2[i][j];
			cnt1[i][j]=cnt2[i][j]=0;
		}
		solve(lt,rt,L,mid),solve(lt,rt,mid+1,R);
	}
}
int main(){
	freopen("star.in","r",stdin);
	freopen("star.out","w",stdout);
	read(n),read(m);
	for(int i=1;i<=n;++i)
	for(int j=1;j<=m;++j){
		s[i][j]=gc();
		while(s[i][j]<'a'||s[i][j]>'z') s[i][j]=gc();
	}
	for(int i=1;i<=n;++i)
	for(int j=1;j<=m;++j){
		l[i][j]=(s[i][j-1]==s[i][j]?l[i][j-1]:j);
		u[i][j]=(s[i-1][j]==s[i][j]?u[i-1][j]:i);
	}
	for(int i=n;i;--i)
	for(int j=m;j;--j){
		r[i][j]=(s[i][j+1]==s[i][j]?r[i][j+1]:j);
		d[i][j]=(s[i+1][j]==s[i][j]?d[i+1][j]:i);
	}
	solve(1,n,1,m);
	cout<<ans;
	return 0;
}
